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Abstract. We consider the family of two-sided initial conditions for TASEP which, as the left 
and right densities (p_,p+) are varied, give rise to shock waves and rarefaction fans - the two 
phenomena which are typical to TASEP. We provide a proof of Conjecture 7.1 of [2U] which 
characterizes the order of and scaling functions for the fluctuations of the height function of 
two-sided TASEP in terms of the two densities p_ , p+ and the speed y around which the height 
is observed. 

In proving this theorem for TASEP we also prove a fluctuation theorem for a class of corner 
growth processes with external sources, or equivalently for the last passage time in a directed 
last passage percolation model with two-sided boundary conditions: p_ and I — p+. We provide 
a complete characterization of the order of and the scaling functions for the fluctuations of this 
model's last passage time L(N,M) as a function of three parameters: the two boundary/source 
rates p_ and 1 — p+, and the scaling ratio -y 2 — M/N. The proof of this theorem draws on 
the results of [12] and extensively on the work of [2] on finite rank perturbations of Wishart 
ensembles in random matrix theory. 



1. Introduction and Results 

We study the fluctuations of height function for the Totally Asymmetric Simple Exclusion 
Process (TASEP) - a stochastic process of great interest due to its wide applicability and mathe- 
matical accessibility. Under hydrodynamic scaling this height function is the integrated solution 
to the deterministic Burgers equation [15]. This hydrodynamic limit is sensitive to the initial 
conditions of TASEP. It is of great interest to determine how the initial condition of TASEP 
affect the random fluctuations of the height function. Ultimately one would like to have a dic- 
tionary between initial conditions of TASEP and the resulting orders of the fluctuations of the 
height function, along with the scaling functions and correlation structures. This paper serves 
to lay some groundwork for understanding the phenomena which figure into this dictionary. The 
two phenomena which must be considered in TASEP are shocks and rarefaction fans. We study 
the simplest family of initial conditions which give rise to both of these phenomena. These ini- 
tial conditions are simply Bernoulli independent at each site x, with density p_ for x < and 
p+ for x > 0. We study the fluctuations of the height function or equivalently the current for 
these two-sided initial conditions. We solve an important conjecture of Prahofer-Spohn [2U] (see 
also |12j). Understanding the fluctuation theory for two-sided TASEP provides the logical link 
between the well developed theory for equilibrium initial conditions (p_ = /»+) [12] and step 
initial conditions (p_ = l,p+ = 0) [14j . Two-sided TASEP interpolates between systems which 
are in equilibrium and systems which are entirely out of equilibrium. Our analysis shows how 
this interpolation occurs. The main result, Theorem ll.l[ was first conjectured in |20j based on 
a scaling theory and analogous results for the PNG model [1]. Figure Q] illustrates the results of 
this paper - it shows how the order of and scaling functions for the fluctuations of the height 
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function for TASEP depend on the observation location with respect to shocks and rarefaction 
fans. 

The proof of our main results makes use of the aforementioned result of [12] for the critical 
point (equilibrium /?_ = p + = p and y = 1 — 2p). For every other set of initial conditions 
the proof relies on the main result of [2], a paper about the largest eigenvalue of finite rank 
perturbations of Wishart (sample covariance) random matrix ensembles. The connection between 
these two, seemingly disparate mathematical construction (TASEP and Wishart ensembles) is 
due to |14j and is facilitated through an intermediate random process known as directed last 
passage percolation (LPP). Our results about fluctuations of currents (or height functions) for 
TASEP result from equivalent results for fluctuations of the last passage time for a directed last 
passage percolation model with two-sided boundary conditions (given in Theorem [L3]) . 

Returning to the model, TASEP is a Markov process rjt with state space rjt € {0, 1} Z . For a 
given t £ M + (time) and x € Z (site location), we say that site x is occupied at time t if rjt(x) = 1 
and it is empty if rjt(x) = 0. Given an initial configuration j]q of particles, the TASEP evolves 
in continuous time as follows: Each particle waits independent exponentially distributed times 
and then attempts to jump one site to its right; if there already exists another particle in the 
destination site, the particle does not move and its waiting time resets (see |15t [TB] for a rigorous 
construction of this process). In equilibrium (or stationary) initial conditions (parametrized by a 
number p £ [0, 1]) the t]q{x) are independent Bernoulli random variables with P{t)q{x) = 1) = p. 
In step initial conditions t]q(x) = 1 for all x < and zero otherwise. Finally, in two-sided initial 
conditions (parametrized by a left density p_ and a right density p + ) T]q(x) are independent 
Bernoulli random variables with P(t]q(x) = 1) = p_ for x < and P(j]q(x) = 1) = p+ for x > 0. 

A natural and important quantity to study in TASEP is the current of particles past an 
observer moving with speed y. It is defined as J y t t t = number of particles to the left of the origin 
at time zero and to the right of yt at time t minus number of particles to the right of the origin 
at time zero and to the left of yt at time t. The current encodes the same information as the 
height function ht(j) (which we will define in equation ([6])): 

ht(j)_-j 
2 

For equilibrium TASEP with density p, the law of large numbers and central limit theorem [S] 
states that 

lim — = (p(l — p) — yp) , almost surely, (2) 

t^oo I 

lim JyM ~ 5 JyM) = N(0, Dj), (3) 

where iV(0, Dj) is a normal with variance 

Dj = p(l-p)\(l-2p)-y\. (4) 

For every velocity aside from y = l—2p, current (and height function) fluctuations are Gaussian 
of order i 1 / 2 . However for a single critical velocity the central limit theorem of [9] is degenerate 
as the fluctuations are of a lower order than t 1 / 2 . In terms of the hydrodynamic limit, this 
velocity corresponds to the slope of the characteristic line for Burgers equation. Heuristically 
this is the speed at which the initial condition fluctuations travel. Therefore, at any other speed, 
the current will depend on more initial conditions than just that localized to the origin - it is 
this that ensures the t 1 / 2 fluctuations and Gaussian scaling function for other velocities. At the 
critical speed, the initial environment's fluctuations are of constant order, and only the dynamic 
fluctuations (those due to the actual TASEP process) are felt. These dynamic fluctuations are 
of central importance to understanding KPZ universality. At the critical speed y = 1 — 2p the 
fluctuations are of order t 1 / 3 and converge, under suitable centering and scaling to a distribution 
function related to the Tracy Widom GUE distribution |2U[ I12j. Rewriting expression (1.14) of 
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Figure 1 . Depiction of three types of TASEP (particles move right) time evolu- 
tion and identification of different regions of fluctuations for corresponding height 
functions. The top diagram depicts the phenomena of a rarefaction fan. The 
height function fluctuations seen by an observer moving at a speed so as to be: 
outside of the fan (outside the dashed lines) will be of order i 1 / 2 and Gaussian 
(denoted Gx)', inside of the fan (inside the dashed lines) will be of order i 1 / 3 and 
Tracy- Widom GUE, (denoted Fq); on the edge of the fan (on the dashed lines) 
will be of order i 1 / 3 and Tracy- Widom GOE 2 , (denoted Fx). Likewise, the middle 
diagram depicts the phenomena of a moving shock. The height function fluctu- 
ations seen by an observer moving at a speed so as to be: on the shock (on the 
dashed line) will be of order t 1 / 2 and "Gaussian squared", (denoted (Gi) 2 ); away 
from the shock (off the dashed line) will be of order t 1 / 2 and Gaussian, (denoted 
Gx)- The bottom diagram depicts equilibrium initial conditions. The height func- 
tion fluctuations seen by an observer moving: at the critical speed y c = l — 2p (on 
the dashed line) will be of order i 1//3 , (denoted Fn and corresponding to what 
[T2] call Fq); at all other speeds (off the dashed line) will be of order t 1 / 2 and 
Gaussian, (denoted Gx)- 

[12] in terms of the current J y t,t, with y = 1 — 2p, their w = 0, and x = p(l — p)i the result shows 
that 

£™ P (~p73^7? ^ x ) = Fi,i{x;0;0), (5) 

where Fx t x(x;0;0) = j^(Fo(s)g(x,0)). The g(x, 0) is a scaling function given in their equation 
(1.18). See Section [1.21 for an overview of how our notation translates into the notation used in 

mini]. 



In TASEP started with step initial conditions there are no fluctuations in the initial envi- 
ronment and consequently for every velocity y S (—1,1) the current has fluctuations from the 
dynamics of order t 1 / 3 and with scaling function which corresponds to the Tracy Widom GU E 
distribution [H] (which we write as Fq so as to be in line with the notation of [2]). In terms of 
the hydrodynamic limit, the range of speeds y G (—1,1) corresponds to the entire rarefaction fan, 
and on the fan, the initial environment fluctuations are not felt versus the dynamic fluctuations. 

Drawing on the heuristics about the fluctuations along flat and fanned regions in the hydro- 
dynamic limit, as well as based on a scaling theory and previous work of [I] for the PNG model, 
Prahofer and Spohn [20] conjectured that these two fluctuation theorems (for equilibrium and 
step initial conditions) arise as cases of a complete fluctuation theory for two-sided TASEP (see 
Figured]). In their Conjecture 7.1, Prahofer and Spohn claimed that the critical point in [9] of 
i 1 / 3 fluctuations for equilibrium TASEP becomes a critical window (representing the region of 
the rarefaction fan) as p_ is increased and p + decreased. Ultimately, as p_ = 1 and p + = 
the critical window of velocities equals the interval (—1,1) as showed in [II]. Likewise they 
conjectured Gaussian behavior outside of this window, as well as in the case where p_ < p + . 

Previous to this paper, part of the Prahofer-Spohn conjecture had been proved via random 
matrix techniques in both the papers of Nagao and Sasamoto [IB] and Baik, Ben Arous and Peche 
[2]. Both papers essentially dealt with the case of p+ = and any p_ £ [0, 1]. Our results are 
dependent on coupling arguments which allow us to bootstrap these boundary cases into every 
type of two-sided initial condition except that critical one (which is dealt with via the result of 

D2I)- 

Much effort has been devoted to understanding the analogous picture for ASEP, where particles 
may move to either the left of the right but are still subject to the exclusion rule. Progress in 
this direction was made in |11| [8j [9j [TO] in the early 1990s. Very recently, due to the efforts of 
Tracy and Widom [291 EH (SB ETJ ) , Derrida and Gerschenfeld [7J, Balazs and Seppalainen [6l[5], 
Quastel and Valko [22j . Mountford and Guiol [17] significant progress has been made in answering 
this question in the general ASEP. The work of Baik and Rains [4|, Prahofer and Spohn |21j . 
and Imamura and Sasamoto [13] in the context of the closely related PNG model has also been 
important in understanding the theory of fluctuations. 



1.1. Results. The main result of this paper is a complete proof of [20J Conjecture 7.1 - our 
Theorem 11.11 

Following [20], assign to a TASEP configuration rjt(j) the height function 



htU) 



'2Jo,t + £ii(l-2»fc(i)), j>l, 
2J ,t, j = 0, 

,2Jo, t -Et,+i(l-2r/tW), j<-l, 



(6) 



Recall that J^t is defined as the number of particles which have crossed the bond (0, 1) up to 
time t. For \y\ < 1 denote 



lim -h t {[yt]) = h(y), 



(7) 



which exists almost surely due to the law of large numbers established via the hydrodynamic 
theory [22 [23] • This limit h depends not just on y, but also on p_ and p + as follows: 
If p- < p + , then 



(l-2p_)i/ + 2p_(l-p_) for y<y 
(1 - 2p + )y + 2p + (l - p + ) for y > y, 



(8) 



with y c = (p+(l - p+) - p_(l-p_))/(p + -p_) = 1 - (p_ +p+). 



If p- > p + , then 

r (l-2p_)y + 2p_(l-p_) for y<l-2p_, 

^(y) = U(y 2 + l) for l-2p_ <y< l-2p+, (9) 

k (l-2p + )y + 2p + (l-p + ) for 1 - 2p+ < y. 

We prove the following (note that in parenthesis we record the distribution names as used in 
For an illustration of the results below see Figure [TJ 



Theorem 1.1 ([20] Conjecture 7.1). (Fq) Let either p^ < p + , y > y c , and —p + < y < 1 — p + , 
or p_ > p + , y > 1 — 2p + , and —p + < y < 1 — p + . Then 

lim P p _, p+ (th(y) - h t ([yt}) < (4p+(l - p+)(y - 1 + 2 P+ )t) 1 ' 2 x) = G x (x). (10) 

t — >oo \ / 

Let either p_ < p + , y < y c , and —p- < y < 1 — p-, or p_ > p + , y < 1 — 2p_, and 
—p- < y < 1 — p-. Then 

lim P p _, p+ (th(y) - h t ([yt}) < (4p_(l - p_)(-y + 1 - 2p^)t) 1 ' 2 x) = G x {x). (11) 



(Fq) Let p- < p + and y = y c , then 

m P p _, p+ (th(y) - h t ([yt}) < ((p + - p-)*) 1 /^) = G 1 ((4p + (l-p + ))- 1 / 2 x)G 1 ((4p_(l-p_))- 1 / 2 : 

(12) 

(Fgue) Let p- > p + and 1 — 2p_ < y < 1 — 2p + . T/ien 

lim P p _, p+ (tfc(y) - ^([yt]) < 2- 1 / 3 (l - y^H^x) = F (x). (13) 
(Fq OE ) Let p- > p+ and either y = 1 — 2p- or y = 1 — 2/? + . T/ien 

lim P p _, p+ (^(y) - ^([yt]) < 2- x / 3 (l - y^fh^x) = F 1 (x;0). (14) 

t — >oo \ / 

(Fq) Let p- = p = p+ and y = 1 — 2p. Then 

lim P p _, p+ ft%)-h t (H) < 2- 1 / 3 (l -yifWx) =F 1 , 1 (x;0;0). (15) 



t^oo 



Remark 1.2. There is one difference between what we prove in Theorem 11,11 and what is stated 
in Conjecture 7.1 of [20] which is that for the case of the Gaussian scaling limit, by virtue of the 
fact that our proof goes by way of a mapping with directed last passage percolation, there are 
certain parts of the Gaussian region (with respect to y,p_,p+) for which our methods do not 
work. 

In this model random weights Wij are associated to each site (i,j) in the upper-right corner of 
Z 2 (with i > and j > 0). The weights are usually independent, and often exponential random 
variables, or geometric random variables. Every directed (up/right only) path ir from (0,0) to 
(N,M) then has weight T(ir) = Yl(i,j)&r w i,ji the sum of all weights along the path. The last 
passage time from (0, 0) to (TV, M) is the maximum path weight over all directed paths: 

L(N,M)= max T(vr), (16) 

n:(0,0)^(N,M) 

The statistics of L(N,M) are dependent on the choice of distribution for the random weights, 
and in certain cases related to eigenvalue statistics for random matrices. 

The current fluctuations for TASEP with step initial conditions were determined by identifying 
the height function with a corner growth model whose growth times correspond with the last 
passage times for a specific LPP model with independent rate one exponential random weights 
Wi j for i, j > and boundary weights Wij = 0, for i = or j = [14] . Theorem 1.6 of [TJ] shows 



that as M, N — > oo such that M/N is in a compact subset of (0, oo), L(N, M) (as defined by the 
above weights) is approximated in distribution by 

where xo is distributed as a Tracy Widom GU E distribution. The first term gives the asymptotic 
average for L(N, M) and the second term shows that the fluctuations scale like M 1//3 and have a 
well understood scaling function. Via the height functions mapping, these results translate back 
into the current fluctuations for TASEP with step initial conditions. The corner growth model 
height function is exactly the random interface bounding the growth region. 

This theorem was proved by using the tools of generalized permutations, the RSK correspon- 
dence and Young Tableaux, to relate the distribution of the last passage time to the distribution 
of the largest eigenvalue of a Wishart ensemble, whose statistics are known to follow the Tracy 
Widom GUE distribution [T3] . 

By analogy our method of proof is to first relate two-sided TASEP to a LPP model, which we 
appropriately call LPP with two-sided boundary conditions (see equation (I18D for a definition), 
and then to relate the statistics of the last passage time for that model to the statistics of 
eigenvalues of already studied random matrices. The first mapping is already found in |20| and 
relies on Burke's Theorem (we review this mapping in Section [3]). LPP with two-sided boundary 
conditions is not directly connected to a random matrix ensemble, however we can realize its 
last passage time as the maximum of last passage time for a pair of coupled LPP with one-sided 
boundary conditions (see (I26D for a definition). The last passage time in such a one-sided LPP 
models is related (see Section 6 of [2|) to the largest eigenvalue of Wishart matrices with finite 
rank perturbations. In fact, the phase transitions, with respect to the magnitude of the finite 
perturbation, which are discussed in [2] correspond exactly to the transitions between different 
orders of and scaling functions for the height function of two-sided TASEP. By a set of coupling 
arguments, and using the results of [2] and [12] we provide a proof of Theorem 11.11 In proving 
Theorem II .11 we reprove the fluctuation results for step initial conditions as well as for equilibrium 
initial conditions (except at the critical point). We also show that these two results arise from a 
much more complete picture (see Figure [T] for an illustration of this). 

As noted above, the proof of this theorem relies on understand the fluctuations of the last 
passage time in a LPP model with two-sided boundary conditions. The specific LPP with two- 
sided boundary conditions which we will devote much of this paper to studying has three different 
types of independent exponential weights Wij: 



w 



exponential of rate 7r if i > 0,j = 

exponential of rate r\ if % = 0, j > 

exponential of rate 1 if i > 0, j > 

zero if i = 0,j = 



(18) 



In the later part of Section [2] we will allow for more general boundary condition where a finite 
number of columns and rows can have different (though uniform within the column or row) rates. 
The LPP with one-sided boundary conditions is defined similarly using the weights in equation 
(I26p . At this point it is worth remarking that changing the distribution of a finite number of 
weights does not have any affect on the asymptotic fluctuations of the last passage time (see 
Lemma 13. ip . 

The statistics considered in this paper are the last passage times L2(N, M) (we use a subscript 
2 to denote two-sided), from (0,0) to (N,M), as N and M go to infinity together such that 
M/N = 7 2 . Note that N denotes the number of columns and M the number of rows. Such 
statistics can be parametrized in terms of the two boundary condition rates tt and n, as well as 
the scaling parameter 7. It is worth keeping in mind that the boundary rates 7r and rj correspond 



with the TASEP densities p_ and 1 — p + , and the scaling parameter 7 corresponds (in a slightly 
more complicated way) with the TASEP velocity y. 

With this connection in mind, we completely characterize both the order and the scaling 
functions for the fluctuations of the last passage time of LPP with two-sided boundary conditions 
in terms of the three parameters tt, 77, and 7. As noted before the main result we appeal to in 
this paper is from [2] (extended to the case 7 < 1 in [19]) which classifies the fluctuations of the 
largest eigenvalue of complex Wishart ensembles with finite rank perturbations. There is a single 
critical point which does not yield to our method of argument, but this corresponds exactly with 
the critical point considered in [12]. Using these two results and coupling arguments we prove 
our LPP with two-sided boundary conditions classification theorem: 

Theorem 1.3. 

(1) For 7 G (0, oo), and M/N —* j 2 , then for tt, 77 such that tt > j^— and 77 > (the GUE 
region ) 



P(L 2 (N, M) < (1 + 7 - x ) 2 M + ( 1+ ^ 4/3 ' M l ^x) -> F (x), 



(19) 



where Fq(x) is the Tracy Widom GUE distribution function. 

(2) For 7 G (0, 00), and M/N — > 7 2 ; then for tt, 77 such that tt > and r\ = or tt = 

and n > (the GOE 2 region), 

P(L 2 (N, M) < (1 + 7 ~ 1 ) 2 M + (1+ 7)4/3 M^ 3 x) -> F^x), (20) 

7 

where Fi(x) is the square of the Tracy Widom GOE distribution function. 

(3) For 7 G (0, 00), and M/N -> 7 2 , then for tt = 1/(1 +7) and 77 = 7/(1 + 7), 

P(L 2 (N,M) < (1 + 7 ~ 1 ) 2 M+ (1 + 7) ^ M 1//3 x) -> Fi,i(x;0;0), (21) 

7 

where Fi i(x;0;0) is the same distribution as what |12] refer to as Fq(x). 

(4) For 7 G (0,oo), and M/N — ► 7 2 , i/ien /or 7r, 77 such that tt < 1/(1 + 7) and rj > 

(the G (tt controlled) region), 

(/ — 1 2 \ / —2 2 \ 1/2 \ 

L 2 {N,M) < (tt-i + JLJLj N+ [tt- 2 - prj^ j ^ 1/2 -J ^ (22) 

where G\ = erf . 

Likewise for tt,7] such that 77 < 7/(1 + 7) and 77 < 7r ^ 1 _^i r 2 ^ + ^_ 2 ('t/ie G (7/ controlled) 
region), 

p [l 2 {n, m) < (V + m + (V - t^Stiji) V2 Ml/2 ^) - (23) 

where G\{x) = erf(x) 

(5) For 7 G (0, oo), and M/N — ► 7 2 , i/ien /or 7r, 77 swc/i i/iai vr + 77 < 1 and 77 = 7r ^ 1 _^Z r 2 ^_ ) _^_ 2 
fi/ie G 2 Zmej, 
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Figure 2. Fluctuation diagram for 7 = 1. Note that all G\ and (G1) 2 regions 
have M 1 / 2 order fluctuations while all other regions have M 1 / 3 order fluctuations. 
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(24) 



)• (25) 



It is worth noting that there are many ways to write the expressions above, and our choices 
are to facilitate the greatest ease in our proofs. 

This type of fluctuation classification picture has been previously discussed in [24j and [3]. 
In fact, in [4j Baik and Rains provide a proof of an analogous fluctuation classification result 
for a closely related particle system model known as the polynuclear growth model. Recently 
[3] studied two-speed (different though related to two-sided) TASEP and proved a fluctuation 
classification theorem for that model. 

1.2. Notation. In a paper such as this which connects two different lines of thought, it is easy 
to become lost in the disparity between notations. We will adopt notation in the style of [2] 
throughout, and when making connections with distributions as found in papers such as [20| . 
|12j . [3J, we will take care to make note of the alternative notation used in those contexts. In this 
Section we define all of the distributions which we will encounter herein and provide references 
for their previous use and definition. 

(1) Grfc(x) is a family of distributions defined in [2] Definition 1.2 and Lemma 1.1. It represents 
the distribution of the largest eigenvalue of a k x k GUE. From this representation it is 
clear that G\{x) = erf(x), the standard Gaussian distribution function. 
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Figure 3. Fluctuation diagram for 7 = 2. Compared with Figure [TT| the effect 
of changing 7 is that the region of M 1 / 3 fluctuations has shifted up and to the left 
along the anti-diagonal. 

(2) Fj(x; xi, . . . , xj) is a family of distributions defined in [2] Definition 1.3. In the case when 
the Xj = for all j, these distributions coincide with those from ||2j Definition 1.1. Of note 
is Fq(x) which is often written as Fgue, the GUE Tracy Widom distribution function, 
and Fi(x; 0) which is often written as Fgoe{x) 2 , where Fgoe is the GOE Tracy Widom 
distribution function. 

(3) Fjj(x; xi, . . . , xj; y±, . . . , yj) is a family of distributions which we conjecture come up 
in LPP with thick two-sided boundary conditions. The only member of this family for 
which we know the correct definition is Fi t i(x; 0; 0) which corresponds to the distribution 
denoted by Fq in [12j . As of yet, we do not know how the other distributions should be 
defined. 

1.3. Outline. The main theorems (Theorem 11.11 and Theorem 1 1 . 3[) have already been recorded 
above in this section. Section [2] provides an intuitive sketch of the proof for Theorem 11.31 Section 
[3] explains the connection between the LPP with two-sided boundary conditions and the TASEP 
with two-sided initial conditions as well as briefly sketches how to translate the result of Theorem 
ll.3l into a proof of Theorem ll.il The final Section 0] gives the full proof of the two main theorems, 
complete with the necessary technical lemmas for the coupling arguments. 
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Figure 4. Fluctuation diagram for 7 = .5. Compared with Figure [TTTl the effect 
of changing 7 is that the region of M 1 / 3 fluctuations has shifted down and to the 
right along the anti-diagonal. 
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2. Fluctuations in Last Passage Percolation with Boundary Conditions 

We start this section by reviewing the result of [2] which relates directed last passage percola- 
tion with boundary conditions to finite rank perturbations of Wishart ensembles. We then apply 
these results to prove Theorem 11.31 which fully characterizes the fluctuations of last passage times 
in terms of boundary conditions and the ratio M/N = 7 2 . Using coupling arguments, supple- 
mented in one case by the result of [12] we provide both the order and the scaling function for 
these fluctuations. Using the exact same arguments but fully taking advantage of the scope of 
the results of [2] , we prove almost all of the cases in Partial Theorem 12.11 In this section we will 
only sketch our proofs, which can be found in entirety in Section [H 

2.1. LPP with one-sided boundary conditions. Consider a directed last passage percolation 
model with one-sided boundary conditions defined as follows: 



exponential of rate 77 
exponential of rate 1 



Hi = oj > 

if i > 0, j > (26) 
zero if i ; > 0,j = 0. 

Let Li(N,M) denote the last passage time from (0,0) to (N,M) (for this LPP model with 
one-sided boundary conditions, but also for any LPP model with thick one-sided boundary con- 
ditions). Then the distribution of Li(N, M) is related to the distribution of the largest eigenvalue 
of the normalized covariance matrix j^XX' where X is N x M and each column is drawn (inde- 
pendent of other columns) from a complex ^-dimensional Gaussian distribution with covariance 



matrix E. The matrix £ has eigenvalues all equal to one aside from a single one, which is Zi = 77 . 
Depending on the value of r/ _1 , L\(N, M) behaves differently. 

The following theorem is adapted from Theorem 1.1 of [2] and the extension to all 7 E (0, 00) 
given in [19], as applied to the one-sided boundary condition LPP. The connection between the 
largest eigenvalue and the LPP with one-sided boundary conditions give above is explained in 
Section 6 of [2] and is briefly rehashed in Remark 12.21 



Proposition 2.1. With L±(N,M) defined as above, as M,N — ► 00 while M/N = 7 2 is in a 

compact subset of (0, 00), the following hold for any real x in a compact set. 

(1) Whenrj> 

P \Lt (N, M) < (1 + 7 - x ) 2 M + ^ + ^ )4/3 M 1 / 3 ^ ^Fo(x); (27) 

(2) W7»enr/ = I ^, 

P ^L X (JV, M) < (1 + l^ fM + ^tl^ M 1 / 3 ^ ^(x); (28) 

(3) Whenr ] <-^, 

P ( Ll {N, M) < (V + ^L^j M + (r?- 2 - J^"^ ) ^ M 1 / 2 ^ - (29) 

Remark 2.2. The connection between last passage time in LPP with one-sided boundary con- 
ditions and the largest eigenvalue of the spiked Wishart ensemble was observed in [5]. The 
connection is not via an exact map but rather a result of the equality of asymptotic distribution. 
Proposition 6.1 of [2] records this fact and explains how a modification of the argument in [T5] 
can be used to prove this. 

An intuitive explanation for the cutoff of n~ l = 1 + 7 _1 in terms of a simple calculus problem 
of maximizing the law of large numbers for LPP paths forced to travel a specific fraction of the 
way along the left column can be found in Section 6 of [2]. 

2.2. LPP with two-sided boundary conditions. Presently we turn our attention to the LPP 
with two-sided boundary conditions as define in equation (|18j) : On the left-most column there 
are exponential weights of rate 77, on the bottom-most row there are exponential weights of rate 
7r at the origin there is a weight of zero, and for all strictly positive lattice points the weight is of 
rate one. Define respectively X(N,M) and Y(N,M) as the coupled last passage times of paths 
which have taken the first step to the right and the first step up (respectively). One should be 
careful to note that we are not conditioning on the location of the optimal path, but rather, for 
each configuration of weights, defining X to be the length of the optimal path which first goes 
right, and Y the length of the optimal path which first goes up. It is clear then that X and Y 
are coupled, dependent and that 

L 2 (N, M) = max (X(N, M),Y(N, M)) . (30) 

Consider now the marginals of X and Y and observe that each of these marginals is of the type 
of the last passage time for a LPP model with one-sided boundary conditions. The boundary 
conditions for Y are exactly as above (rj weights and an N by M region). However, for X, the 
boundary conditions are it weights and an M by JV region (note that the region has been flipped 
in order to conform with the setup for Proposition I2.ip . From this observation we can apply 
Proposition 12.11 to completely characterize the marginals of the joint distribution for the pair 
(X,Y). Note that while X and Y are not exactly of the form of a last passage time for a LPP 



with one-sided boundary conditions, they only differ by a finite number of weights and therefore 
have the exact same asymptotic statistics via Lemma 13. 11 



Proposition 2.3. With X(N,M) defined as above, as M,N -> oo while M/N 
compact subset of (0, oo), the following hold for any real x in a compact set. 

(1) WhemOj^, 



P \X(N, M) < (1 + 7) 2 iV + 

(2) Whemr = j^, 

P (x(N, M) < (1 + 7) 2 iV + 

(3) Whenir <T ^ ; 

P lx(N,M) < (V x + 7T _ 1 1 ^ ]n + 



(l+ 7 -i)V3 



^l +7 -l)4/3 



F n (x) 



7 is in a 



(31) 



Fi(x) 



(32) 



7T- 2 7 2 

(vr- 1 - l) 2 



1/2 



GUx). 



(33) 



Proposition 2.4. With Y(N, M) defined as above, as M, N — ► oo m/iife M/N = 7 2 is in a 
compact subset of (0, oo), i/ie following hold for any real x in a compact set. 



1) MenpJ-, 



P I y(JV, M) < (1 + 7" 1 ) 2 M + 



7 



(2) Whenrj 



P I y(JV, M) < (1 + 7 _1 ) 2 M + 



;i+7) 4/3 

7 



Fi(x); 



(34) 



(35) 



(3) fci,<^ 



p y(iv, m) < 



-1 —2 

-1.7/ 7 
r? + 



V 



-i 



M + ?? 



r] 2 7 2 
(r?- 1 - l) 2 



1/2 



M 1//2 x ] ^Gi(x); (36) 



We assume that both n and ir are both between zero and one. In fact, it is clear from our proofs 
that the order and fluctuation of our two-sided last passage time L2 (N, M) for parameters rj and 
7T is that same as that for parameters rj A 1, it A 1. Thus it suffices to consider only n, ir £ [0, l] 2 . 

For each of X and Y there are two regions of different fluctuation orders, and one critical 
point which has 1/3 order fluctuations. We call the point (vr,?]) = (1/(1 +7), 7/(1 +7)) the 
critical point for the pair tt, ry. If vr < then the X fluctuations are of order 1/2, and likewise 
if rj < y^— then the Y fluctuations are of order 1/2, whereas in the complementary cases, the 
fluctuations are of order 1/3. By comparing leading (law of large number) terms in Propositions 
12.31 and 12.41 we see that if either of these two inequalities hold, then the fluctuations must be of 
order 1/2. In this case, either the leading term for X or Y clearly wins, in which case ^(iV, M) 
has the leading order behavior and Gaussian fluctuations of the winner random variable, or the 
two random variables have the same leading terms. The second case, or equal leading terms, 
occurs when 



TT H ; 



7T 



N 



-1 v S 2 



v 



-1 



M. 



(37) 



V= —i , .,_ 2 - (40) 



In this case the fluctuations will remain of order 1/2, but will behave as the fluctuations of two 
independent normal random variables (what we call G 2 ). 

If 7 = 1 there are two solutions to equation ([37|) . One is rj = tt and the other is rj = 1 — w. 
Since we are only considering the Gaussian region, the anti-diagonal solution is of no interest, 
and we find that we have G 2 density for our fluctuations if r\ = tt and rj < 1/2. 

For 7^1, the solution set is a little harder. Recall M = j 2 N and using this we can factor 
out N from both sides giving 

Applying the change of variable tt — ► 1 — tt we find that it suffices to solve 

^-' + 5F^t) = (^- , + ^t). < 39 > 

and change the solution back to our original variables. 

This again has the solution rj = it. Solving for the other solution and then changing variables 
back we get 

TT 

7r(l - 7~ 2 ) + 7" 

By plugging in the critical point (1/(1 + 7), 7/(1 + 7)) it is easy to see that this G 2 curve is 
continuous between the origin and the critical point, though only linear for 7 = 1. 

We have, so far, only accounted for the regions where tt < or r\ < j+z- There are four other 
regions to consider which correspond to replacing the or with an and, and the less than sign with 
either equality, or a greater than sign. In each of these cases the leading term is independent 
of tt and rj and equals (1 + 7 _1 ) 2 M. The fluctuations of X and Y are both or order 1/3, so 
those of L2(N,M) are as well. To determine the scaling functions finer coupling arguments are 
necessary. For instance when tt > and 77 > jj^—, the last passage path for X and for Y can 

be compared to the analogous random variables X and Y, for last passage paths for a coupled 
LPP model with two-sided boundary conditions all of rate one (i.e. 7r, 77 = 1) such that pointwise 
X > X and likewise Y < Y. The maximum of X and Y equals L(N,M), the last passage time, 
and we show that X and X, and likewise for Y and Y, under appropriate centering and scaling 
converge to the same distribution respectively. This implies, via our Lemma 14.21 that the scaled 
and centered random variables in fact converge in probability and hence that their maximums 
converge in probability. This means that the maximum of X and Y behaves just like a regular 
last passage time which is known to have GUE scaling function. A similar coupling shows that 
the scaling when either tt = and r\ > j^— , or tt > j^- and rj = j]-^ the scaling function 
behaves like that of a single critical last passage time for a LPP model with one-sided boundary 
conditions. 

Determining the scaling function at the critical point is a harder problem. One may identify it 
as the maximum of two F\ distributions, coupled as X and Y are coupled. This characterization, 
a priori, yields a tight family of random variables. However, how to prove that it converges on 
more than just a subsequence is not immediately clear, and more over it is not clear to what it 
convergences. A posteroi, this characterization is justified since the result of [12] can readily be 
translated into a proof that the scaling function at the critical point is iq^x; 0; 0) (what they 
call F ). 

The results from [2] used in the proof of Theorem 11.31 yield . in fact, a much more general result 
via essentially the same argument. We now define what we call the LPP model with thick two- 
sided boundary conditions in terms of boundary row and column thickness integer parameters 
J, I > 0; two vectors of row and column weight rates tt = (717, . . . , ttj), rj = (rji, . . . , 77/); two 
vectors of row and column convergence rates X = (xi, . . . Y = (yi, . . . ,yi). With these 



(41) 



parameters our model is defined in terms of the following LPP weights (which implicitly depend 
on M and JV): 

' exponential of rate ttj + if i > J, j < J; 

exponential of rate r\i + jrij3 if i < I, j > J] 
exponential of rate 1 if % > I, j > J; 

zero if i < I,j < J. 

To see that this model is a broad generalization of our previously considered two-sided bound- 
ary condition model, take J, I = 1, 7Ti = 7r, rj\ = r) and x\,y% = 0. Corresponding to this 
model we now provide a complete characterization of its asymptotic fluctuations. A number of 
distributions not previously discuss are introduced in this theorem. A full discussion of these 
distributions can be found in Section 11.21 

The coupling arguments given to prove Theorem 1 1 . 3 1 can be easily adopted to this new setting. 
Given the above parameters define, again, two coupled random variables X(N, M) and Y(N, M) 
as follows. X(N,M) is the last passage time from (0,0) to (JV, M) of the set of up/right paths 
which cross through at least one vertex from the set : i = I,j € {0, . . . , J — 1}}. Likewise 

Y(N,M) is the last passage time from (0,0) to (JV, M) of the set of up/right paths which cross 
through at least one vertex from the set ■ i =£ {0, 1}, j = J}. Clearly any up/right 

path from (0, 0) to (JV, M) must go through one and only one of these two regions. Furthermore, 
by virtue of the definition of the last passage time, the maximizing path for X(N,M) and 
Y(N,M) will necessarily go through the points (0, J) and (/, 0) (respectively). Thus we may 
refine our definitions of X(N,M) and Y(N,M) to require passing through these two points. 
Again we see that L2(N, M) = max(X(N,M),Y(N,M)) and just as before [2] provides an 
immediate proof of the following: 

Proposition 2.5. For the vector n fix the set K\ C {1, . . . , J} by 

K 1 = {je{l,...,J}:Tr j = T ^}, (42) 

an define Xk x as the elements of X which correspond to indices in K\ . Further, define tt = 
linXj-g/i ,.. and let k\ be the number of ttj which attain the value tt. 

Then with X(N,M) defined as above, as M,N — > oo while M/N — 7 2 is in a compact subset 
of (0, oo), the following holds for any real x in a compact set. 

(1) fcO^, 

P (^X(N, M) < (1 + j) 2 N + Q+J-p— N^xj -> F (x); (43) 

(2) When 7r = j^, 

P (x(N, M) < (1 + ifN + (1 + ^ 1)4/ V /3 X j ^F lKll (x-X Kl ); (44) 



(3) When 7f< T ^ ) 

P (x(N,M) < (V 1 + ^r^) N+ (V 2 - ( -!r_ 7 i )2 ) ' Nl/2 A ^ (45) 



A similar proposition exists for Y(N, M). Using these two results the same types of coupling 
arguments then apply and give both the orders and the scaling functions for L2 ( JV, M) . As 
before, these coupling arguments break down when both boundary conditions are critical. With 
single width boundary conditions, we appealed to [12J, however in this case no existing argument 



provides a characterization of the behavior in this case. The following partial theorem therefore 
contains a single conjectured equation (f5Tj) whose study seems very difficult. 



Partial Theorem 2.1. 

(1) For vectors tt, r\ fix the sets K\ C {1, . . . , J} and K 2 C {1, . . . , /} by 

K 1 = {j€{l,...,/}:n j = T ^} (46) 

K 2 = {ie{l,...,I}:r H = -?-}. (47) 

1 + 7 

Then define Xk x and Yk 2 as the elements of X and Y which correspond to indices in K\ 
and K 2 respectively. 

For 7 G (0,oo), and M/N — > 7 2 ; then for vectors tt,t] such that ttj > j-jj— for all 
i £ {1, . . . , J} and r\i > for all i 6 {1, . . . , /} then if: 

(a) |#i|=0, \K 2 \=0, 

P(L 2 (N, M) < (1 + 7 _1 ) 2 M + (1+ 7)4/3 M^x) F (x); (48) 

7 

(b) |Ki| >0 ; |K 2 | =0, 

P(L 2 (N,M) < (l+ 7 -l) 2 M + (1+ 7)4/3 M X /3 X ) -> (49) 

7 

(c) |Ki| =0, |tf 2 | >0 5 

P(L 2 (N,M) < (l + 7 - 1 ) 2 M+ (1+ 7)4/3 M^x) ^F|^ 2 |(x;y^ 2 ); (50) 



(d) |Ki| >0, |K 2 | >0 ; 
P(L 2 (N,M) < (l+ 1 - 1 ) 2 M+ {1+l)4/3 M l ^x)^F lKlllK2l (x;X Kl ,Y Kl ] 



(51) 



(2) Define tt = min Jg {i ,...,j}(7Tj)j anc ^ *7 = m i n ie{i,...,/}( r ?«)? ari ^ ^ ^1 ^ e ^ e number of TTj 
which attain the value tt and likewise k 2 be the number of r]i which attain the value fj. 
For 7 € (0, 00), and M/N — ► 7 2 , i/jen /or tt, 77 smc/i i/iaf: 

(a) 7f < 1/(1 + 7) anc ^ ^ > -^ 1 _^i r a - )+ ^_ 2 f/j/ie G (7r controlled) region), 

(/ — \ 2 \ / — 2 2 \ 1/2 \ 

£ 2 (iV,M) < ^T 1 + JV + ^-2 _ ^rj^ j iV 1/2 x J - G fe (x); (52) 

(b) fj < 7/(1 +7) and 77 < - (1 _ 7 _a )+7 -a f^e G C 7 ? controlled) region), 

(/ ~— 1 — 2 \ / ~— 2 —2 \ 1/2 \ 

l 2 (n,m) < (rj- 1 + frnrij M + - (^-1-1)2 ) Ml/2x j ^ ^^) ; ( 53) 

(c) tt + 77 < 1 and fj = ~ {1 _^- i)+ „ r i (the G 2 line), 

P L {N , M) s + -i) , + (fiz£±*g^z«vi) 1/2 (54) 



Gfcifo = = )G k2 {x ). (55) 

V 7 ! - TT + 7 2 7T 7 



Finally let us note two applications of LPP with thick one-sided boundary conditions which 
can be found in pQ. The first application deals with what Baik called traffic of slow start from 
stop in which particles start in the step initial condition of TASEP and have a start-up profile - 
that is, every particle moves slower for its first few jumps, and then returns to jumping at rate 
one. The second application is dual to the first one and is called traffic with a few slow cars in 
which particles always move at a slower rate. In both cases Baik identifies the fluctuation scaling 
limits by using the [2] type results which we have made use of herein. 

In the next section we will give an important application for the LPP with two-sided bound- 
ary conditions model to two-sided TASEP. It is unclear whether the thick two-sided boundary 
conditions model has any similar application to TASEP or related models. 



3. Mapping TASEP to Last Passage Percolation with Boundary Conditions 

In this section we explain the connections between the last passage time in LPP with two-sided 
boundary conditions and the fluctuations of the height function for the two-sided TASEP model. 
Making use of this mapping we explain how the results of Theorem 11.31 imply the results of [9] 
stated in the introduction. Furthermore we briefly explain how this theorem translates into a 
proof of Conjecture 7.1 of [20] (full proof is giving in Section H]). 

We start with a lemma which states that finite perturbations of our LPP model, have no affect 
on the asymptotic behavior of the last passage time. 

Lemma 3.1. Fix some LPP model with weights Wij (independent but not necessarily identically 
distributed) such that 

P ( L{N >^- m < *) - F( X ) (56) 

for M/N — > 7 2 G (0,oo), and for F a non- degenerate probability distribution. Randomly, 
independent of the values of w^j, change a set of these weights to a new set of weights w\ j 
and let L'(N, M) denote the last passage time with respect to the original weight with the newly 
updated weights. Call A the set of changed indices (i,j) and Wa = Yl(i j)£A w i,j + w 'i,j- Then if 
E[Wa] < oo and ifb^ — > oo, we also have 

P (a»< I ). FW , pi) 

Below is an outline of the proof. The full level of details is suppressed since a similar style of 
proof is given for Lemma 14.11 in full detail. 

Proof. Since the total effect of the change of weights corresponding to A has finite expectations, 
the Markov inequality shows that for any e we can find I large enough so that P(Wa > — e - If 
we restrict ourselves to this region of our statespace, then since the bjy goes to infinity, the effect 
of the change of weights is negligible in the limit. Since this is true on all but an e region of the 
state space, we have that the distribution functions are within e of each other in the limit, but 
taking e to zero gives equality. □ 

Recall our definition of the two-sided TASEP model given by the initial conditions of Bernoulli 
with parameter /)_ on the left of zero and with parameter p + on the right. Corresponding to 
the TASEP process started with this random initial condition, we consider the height function 
ht(j) defined in equation ©. Theorem 2.1 of [20] relates the joint distributions for this height 
function to those of the height function for a particular growth model associated with a variant 
on the LPP with two-sided boundary conditions. The weights for this variant LPP are defined 
with respect to two independent geometric random variables Cf and C-. Let Cf be geometric 
with parameter 1 — p + (i.e. P(C+ = n) = p+(l — p+) n ) and C- be geometric with parameter p_ 



(i.e. P(C- = n) = (1 — /3_)p"). The weights are then defined as independent random variables 
with: 

exponential of rate 1 if i,j > 1; 

zero if i = j = 0; 

zero if < z < C+ , KO x 
n>ij = < (58) 

exponential of rate 1 — p+ if i > C+; 

zero if < j < £_ ; 

exponential of rate /?_ if j > 

With respect to these random weights define the last passage time LrN,M). This family of 
random variables is non-decreasing in both N and M. Therefore one can associate to this a 
growth process on the upper-corner and likewise a height process over the number line. Let 
At = {(N,M)\N,M > 1,L(N,M)} be the growth process and let h t be defined so as to satisfy 
At = {(N, M)\2 < N + M < ht(N — M)}. To describe this in words, imagine rotating counter- 
clockwise, the upper corner in which LPP occurs by 7r/4. To each lattice point (labled by j € Z) 
now on the horizontal associate a height ht(j) equal to two times the number of L(N, M) vertically 
above j which are less than or equal to t. Then we have the following: 

Theorem 3.2 (Theorem 2.1 of |20j). In the sense of joint distributions we have 

ht(j) = ht(j) for\j\<h t (j). (59) 

This theorem essentially says that for the height profile which lies above the boundary of the 
rotated upper corner, the two profiles have the same joint distribution. It is worthwhile to recall 
that there is a similar map between the TASEP height function for TASEP with step initial 
conditions and the height function for standard (no boundary condition) LPP p3j. The proof 
of Theorem 13.21 can be found in [20] and essentially amounts to a study of the dynamics of the 
right most particle to the left of the origin, as well as the dynamics of the left most hole to the 
right of the origin. Tagging this particle and this hole we observe that their initial location is 
geometric and using Burke's theorem we find that their waiting times between successive moves 
is exponential with rate relating to the densities p- and p+. Between the tagged particle on the 
left and the tagged hole on the right, particles move according to normal TASEP rules, and hence 
the two height functions evolve with the same dynamics. The dynamics of the boundary of the 
part of the TASEP height function lying in the rotated upper corner is matched by the effect of 
the LPP boundary conditions, and the theorem follows. 

From this theorem we see that the following equality: 

P P -,p+(ht(N -M)>N + M) = P(L(N, M) < t). (60) 

This equality is not of much use to use however because in order to use the results of theorems 
11.31 we must deal with a slightly different LPP model than corresponding to L. However, this 
model and our LPP with two-sided boundary conditions only differ in expectation by a finite 
number of weights (the geometric number of zeros from (_ and £+))• Therefore while it is true 
that there is not exact equality then with P(L2(N, M) < t), from Lemma 13.11 we see that for 
any sort of central limit fluctuation statement with a non-trivial limiting distribution, we have 
equality in the limit. We will abuse notation in the remainder of this section and the next during 
the proof of Conjecture 7.1, and write equality between the TASEP height function probability 
and the probability for the last passage time in the LPP with two-sided boundary conditions 
model. 

It is important to note that Lemma 13.11 only applies if both p_ < 1 and p + > 0. If either of 
these inequalities is violated, then the geometric number of zeros on the boundary will in fact, 
almost surely be infinite. However, in any of these cases, the classification of one-sided LPP then 
readily applies. 



The boundary conditions for L corresponded to having exponentials of rate p- on the left 
boundary and 1 — p + on the right. Therefore in terms of ir and r\ we have tt = 1 — p + and 77 = /?_. 
The critical point for 7r, 77 is -rj— and y^— , therefore we see that the critical point for p~,p+ is 



1+7 1+7' 

P- = P+ = 77— ■ ( 61 ) 
1 + 7 

This corresponds to an equilibrium measure on TASEP with density j^—- 

In the next section we will show how Theorem 11.31 implies an almost complete (all put a few 
regions of the claimed Gaussian region are fully proved) proof of [20J Conjecture 7.1 (our Theorem 
II. ip . From this result we may easily deduce the results of [9] stated in the introduction. 

For simplicity assume r £ [0, 1], as the case r S [— 1, 0] follows similarly. We wish to prove that 

P ( J ^ " Ml ~ P) - rp)t > x)-+ <h{x), (62) 
\ t V^ p{l - p)\{l - 2p) - r\ ~ ) 

where as defined before G\{x) is the standard Gaussian distribution function. 
In the case of p- = p+ = p and r > we can conclude from equation (jlOD that 



P (th(r) - h t ([ry]) < (4p(l - p)(r - 1 + 2p)i) 1 / 2 x) -> G^x), (63) 

where h(r) = (1 — 2p)r + 2p(l — p). Substituting the relationship in equation ([T]) and rearranging 
terms we arrive at the exact result of [9] desired. 

We now briefly explain the approach to proving Conjecture 7.1 from our Theorem 11.31 The 
conjecture deals with height functions. We have provided above the relationship between height 
function distributions and LPP distributions. From equation (|60|) we see that if one is to consider 
P(ht(j) > x) in terms of LPP, you must solve for N = and M = The variables j and 
x both are functions of time t and a speed y. If M/N = has a limit, we call that j 2 . This 
allows us to asymptotically write M (or N) just as a function of time (thus the y dependence 
goes into 7). In the cases we consider we can invert the expression for M in terms of t and get 
an expression for t in terms of M, thus putting us in the form of the limit theorems we proved 
in Theorem 11.31 

4. Proof of Fluctuation Theorems 

In this section we provide a proof of Theorem 11.31 (which easily generalizes to prove Partial 
Theorem 12 . 1 [) and a proof of Theorem 11.11 

4.1. Proof of Theorem 11.31 (Fluctuations of LPP with two-sided boundary conditions). 

The following two technical lemmas provide the basis for the coupling arguments necessary in 
our proof of Theorem 11.31 

Lemma 4.1. If X n > X n and X n D as well as X n D then X n and X n converge in 
probability (i.e. for all e > 0, as n — > 00, P(X n — X n > e) — > 0). Conversely if X n > X n and 
X n => D and X n converges to X n in probability then X n D as well. 

Lemma 4.2. Assume X n > X n and X n =>■ D\ as well as X n ^ D\; and similarly ifY n >Y n 
and Y n =>■ D2 as well as Y n D2. Let Z n = m&x(X n ,Y n ) and Z n = m&x(X n ,Y n ). Then if 
Z n => Z?3, we also have Z n ^ D%. 



Proof of Lemma \4.1\ While it is likely that this lemma is known in the literature, we do not 
know where and hence produce a proof. We prove the first assertion. Fix e > and, from 
the point of contradiction, assume that P(X n — X n > e) > 5 > for an infinite subsequence 
of n's. By restricting to this subsequence and noting that all of the hypothesis of the lemma 
hold under this restriction, we may equivalently assume that P(X n — X n > e) > 5 > for all 



n large. Since X n and X n converge weakly, each sequence of random variables is tight. This 
implies that there exists an M(e) and N{e) such that for all n > N, P(\X n \ > M) < 5/2 and 
likewise P(\X n \ > M) < 5/2. Thus P(\X n \ >Mfl {X n - X n > e}) < 5/2, therefore 

P(\X n \<Mn{X n -X n >e})>5/2. (64) 

Call this event A = \X n \ <Mfl {X n — X n > e}, then conditioned on A, X n > X n + e. For large 
enough n, 

P(X n < t\A) < P(X n < t\A) - P{X n G [t - e,t]\A). (65) 
Now we partition the interval [— M, M] into e size blocks and use the pigeon hole principle to 
show that there exists j E {0, . . . , f^-]} for which 

P(X n G [-M + 6(j - 1), —M + ej]\A) > ± (66) 
Therefore taking t = — M + ej we have 

P(X„G[i-e,i]L4)>^_ (67) 

Therefore 

P(^n < t\A) < P{X n < t\A) - ± (68) 
Multiplying both sides by P(A) and re-writting without conditioning gives 

P(A) 

P(X n <tnA)<P(X n <tnA)-^t. (69) 
That X n y X n also implies that 

P(x n < tn A c ) < P(x n < tnA c ). (70) 

Adding these two inequalities and using the fact that P(A) > 5/2 gives, for all n large enough 

P(X n <t)<P(X n <t)-^. (71) 

This inequality implies, however, that X n and X n can not converge in distribution to the same 
object. This is a contradiction to our hypothesis, so our assumption must be false. That is, 
P(X n — X n > e) must go to zero as n goes to infinity. 

The second assertion is easier. For all e we can find N such that for n> N P{X n — X n > e) < e. 
Set inclusion and partitioning implies that 

P{Xn <t)> P(X n <t)> P(X n <tDX n -X n >e) + P(X n <tDX n -X n <e). (72) 

The first term on the right hand side P(X n < t n X n — X n > e) is less than e by convergence 
in probability. The second term P(X n < t D X n — X n < e) > P(X n < t — e) by inclusion. 
Re-writting the inequality we have 

P{Xn <t)> P{X n <t)>e + P(X n <t-e). (73) 

If t is any continuity point for D, then we can take e to zero and we find that 

D(t) > P{X n <t)> D{t), (74) 

and hence X n weakly converges to D. □ 

Proof of Lemma \4-2[ Applying Lemma 14.11 to both Xn and A^v, as well as Yjy and Y/v we find 
that for any e, large enough N, P(Xn — Xn > e) < e and likewise P(Y/v — Y/v > e ) < e - Therefore 

P(Z n <t)> P(Z n <t)>3e + P(Z n <t-e), (75) 

where the right hand side comes from partitioning the event Z < t with respect to Xn—Xn > e 
and Yjy — Yv > e and noting that three of the partitions have probability trivially bounded by 



e. Taking t to be a continuity point for the limiting distribution for Z (for D3) and taking e to 
zero we get that P(Z < t) = Fr> 3 {t), and hence Z converges in distribution to -D3. □ 

Proof of Theorem \1.3l (Fq) . This result follows from a coupling argument between the X(N, M), 
Y(N, M) variables as well as a second set of last passage times X(N,M), Y(N, M). X(N,M) 
and Y(N, M) are coupled as previously described (they are the last passage times if forced to 
go right (or up) on the first move). Now to define the tilde versions of X(N,M) and Y(N, M), 
for a given realization of weights, divide the boundary weights by their means. This creates a 
new set of weights coupled and pointwise dominated by the original set of weights. For this 
new set of weights define X(N, M) and Y(N, M) as the last passage time if forced right (or up) 
initially. From this pointwise domination of the new weights by the original weights, we see that 
X(N,M) > X(N,M) and Y(N,M) > Y(N, M) pointwise. The advantage of the tilde variables 
is that Z(N,M) = max(X(N,M),Y(N,M)) is the standard (without boundary conditions) last 
passage time. We also know that, asymptotically X(N,M) and X(N,M), as well as Y(N,M) 
and Y(N, M) have the same distribution. We wish to use this information to conclude that 
Z(N, M) and Z(N, M) have the same distribution as well. 

Let us redefine our variables by properly shifting and scaling them, so that they have a limiting 
distribution. Our new X(N,M) is 

X(N,M)-(1 + ^)N 

7 ( 1+7 -l)4/3JVV3 ' 1 > 

and similarly we define X(N, M), Y(N, M),Y(N, M), and Z(N, M) and Z(N, M) is redefined in 
terms of the newly defined variables. Now we have the following setup: X(N, M) > X(N, M) and 
X(N,M) => F as well as X(N,M) => F ; similarly Y(N,M) > Y(N, M) and Y(N,M) => F 
as well as Y(N,M) =>■ Fq. Applying Lemma 14.21 we get that Z(N,M) converges in distribution 
to F . □ 

Proof of Theorem \1.3l (F\). There are two cases which yield to the same argument. Thus we 
prove the case of it > 1/(1 + 7) and 77 = 7/(1 + 7) only. As in the last proof, let X(N,M), 
X(N,M), Y(N,M) denote the suitably centered and scaled random variables, for the sake of 
applying Lemma[0]we define Y(N, M) = Y(N, M). Again we know that X(N, M) > X(N, M) 
and X(N, M) => F as well as X(N, M) => F and clearly the same holds for the Y(N, M) and 
Y(N,M). So by Lemma 14.21 since we know that max(X(A r , M), Y (JV, M)) converges weakly to 
Fx , it follows that max(X (N, M) , Y (N, M)) => F\. □ 

Proof of Theorem (F\^). This follows immediately from the main result of [12]. □ 

Proof of Theorem \1.3\. (G). We prove the first case, when it < 1/(1+7) an d V > ^(i^-^+^-i > as 
the other case has the same proof. In this region of it, r], the leading order term on the expression 
for X(N, M) is larger than that of the expression for Y(N, M). If we renormalize both X(N, M) 
and Y(N, M) by the leading order term for X(N, M) and divide by its fluctuation term, we find 
that X(N, M) converges to a standard normal. On the other hand, since the leading order term 
for X(N, M) exceeds that for Y(N, M), the renormalized Y(N, M) converges to negative infinity. 
This implies that max(X(N, M),Y(N, M)) converges in distribution to a standard normal, just 
like X(N,M). □ 

Proof of Theorem E3 (G 2 ). We couple X(N, M) with a random variable X(N, M) where X(N, M) 
is the last passage time when forced to stay along the bottom edge for at least a specific, deter- 
ministic fraction of the path (we likewise define Y(N, M)). Specifically we define X(N, M) to be 
the last passage time when the path is pinned to the bottom edge until the point 



after which point is forced into the bulk and allowed to follow a last passage path therein. The 
weights accrued along the bottom edge respect a simple central limit (as they are the sums of a 
deterministic number of iid random variables) and the weights accrued after the path is forced into 
the bulk follows the fluctuations theorem for standard exponential last passage times. As these 
two random variables are independent by construction, their means add since their fluctuations 
are of different order (iV 1 / 2 for the bottom and N 1 / 3 for the bulk) the bottom fluctuations win 

~ — 12 

out. Following this idea we find that the mean of X(N, M) is (tt~ 1 + %=j^)N and the fluctuations 
are normal with variance 



To see this, observe that if we center X by the mean and divide by the square root of the 
above variance, we are left with a random variable of the form X(N,M)i + N~ 1 ^X 2 (N, M), 
where X(N,M) converges to a normal, and X 2 (N,M) converges to a GUE. Since the second 
term has a prefactor which goes to zero, we see that P(X(N, M)± + N^ 1 ^ 6 X 2 (N, M) < I) can 
be partitioned into a region of size e where \X 2 (N, M)\ > R and a region of size 1 — e where 
\X2(N, M)\ < R. On the second region, we can replace X 2 (N,M) by R and find asymptotically 
that the probability differs from P(X(N,M) < I) by only e. Taking e to zero gives the desired 
convergence in distribution. 

Therefore if we center X(N,M) and X(N,M) by the same amount and renormalize by the 
same amount we get two random variables which converge to the same distribution, despite the 
first one being almost always larger than the second one. This is one of the pieces we will need 
to apply Lemma 14.21 

We can likewise define Y(N, M) as the last passage time when the path is pinned to the left 
edge until the point 



after which point is forced into the bulk and allowed to follow a last passage path therein. As in 
the prior we see that centered and renormalizing Y(N, M) and Y(N,M) by the same amounts, 
gives to random variables which converge to the same distribution, despite the first one being 
almost always larger than the second one. 

From the relationship between it and r\ we know that the leading order terms for both X(N, M) 
and Y(N, M) coincide. Therefore we can write 




(78) 




(79) 



X(N, M) 
Y(N, M) 



AN + Bir^^X^N, M) + C^N 1 ' 3 X 2 (N M) 
: AN + BrjN^Y^N, M) + C^iV 1 / 3 ^^, M), 



(80) 
(81) 



where 




(82) 



(83) 



(84) 



and C w and C„ are constants 



(which will play no role here). If we consider now Z(N,M) 



max(X(N, M),Y(N, Af )) we find that 




where E\ and E<i are respectively the events 




B *X X {N,M) + C'^N- X I & X 2 {N,M) < x, 



Yi(N, M) + C^N~ 1/6 Y 2 (N, M) < x. 



(85) 



(86) 



As before, because of the iV _1//6 prefactor to the X2(N,M) and i^(iV, M) terms, we can 
condition on these terms being bounded by some large number R, and only cost ourselves e of 
the sample space. Once we have conditioned on these random variables being bounded by R, we 
can conclude that their joint probability is bounded between the product 



P Xi(N,M) < 




C' n N-^ 6 R P fi(N, M) < 




(87) 



P (^Xi{N, M) < yj^-x + C'^N-WrJ P (yi(N, M) < + C^iV" 1 / 6 ^ . (88) 

Taking N to infinity gives P(Xx(N, M) < x)P(Y\(N,M) < x) and taking e to zero, and using 
the central limit theorem to show that X\{N, M) is standard normal, we find that 



P 



Z(N,M) - AN 



< x 




(89) 



Therefore, using the observations made at the beginning of the proof, we can apply Lemma [4. 2 1 
to conclude that the probability distribution of Z(N, M) centered and normalized as above, con- 
verges to the same G 2 distribution (i.e. equation ([89]) holds with Z(N,M) in place of Z(N,M)). 
Working out the coefficients B n and using the relationship between ir and r\ we have our 
desired result. □ 



4.2. Proof of Theorem 11.11 (Conjecture 7.1 from [20]). The following elementary lemma 
will find repeated use in what follows. 

Lemma 4.3. If M = at + bt 1 ' 2 then for large t, 

t = a~ 1 M -a~ 3/2 bM 1/2 + o(M 1/2 ). (90) 
Likewise if M = at + bt 1 / 3 then for large t, 

t = a~ 1 M -a'^bM 1 / 3 + o{M 1 ' 3 ). (91) 

Proof of Theorem \l.ll We provide proofs of only the G 2 and Fq cases of this theorem as the 
G case is analogous to G 2 and the Fx case is analogous to Fq. The Fx t x case of the theorem 
already is proved in [12]. The proves are based on mapping the event that the height exceeds a 
given value to the event that the last passage time is less than that value. Then Theorem 11.31 
applies and gives asymptotic height distribution results. As noted before, the mapping between 
the TASEP height function and the last passage time for our model of LPP with two-sided 
boundary conditions is not exact (as the exact LPP model has geometric numbers of boundary 
zeros) however Lemma 13 . 1 1 ensures that asymptotically all results in our LPP model correspond 
to results for the two-sided TASEP height function. 

G 2 case: Recall that in LPP with two-sided boundary conditions n = 1 — p + and rj = p_. We 
presently assume that n + ir < 1 (p_ < p + ) and y = y c . Recalling that h(y) = (l — 2r])y + 2r](l — n), 
we wish to determine the asymptotic (large t) value of 

P{h t {yt) > ((1 - 2n)y + 2r?(l - n))t - y/(l - vr - r,)t l ' 2 x). (92) 




We will reduce this probability to a probability in the related LPP with two-sided boundary 
conditions, and then use Theorem 11.31 to conclude that this probability is the correct product of 
Gaussian probability functions. 

The first step in translating to a LPP problem is to relate the speed to 7. We may solve for the 
asymptotic value of 7 as a function of the speed y. As shown in Figure 14.21 the height function 
event corresponds to the LPP event where N — M = yt and N + M = ((1 — 2r))y + 277(1 — r)))t — 
\/(l — 7T — rj)t l ' 2 x). From that we find that 

M = 77(1 - y - n)t - ^(l - tt - vjt 1/2 x, (93) 

N = (y - 7777 + 7/(1 - rj))t - - V(l - vr - ^)t l ' 2 x. (94) 

Therefore 

7 2 = lim M/N = ^-y-v) (95) 

From this equation we can solve for y as a function of 7: 

( 7 2 - lMl-ry) 

y= — 27 — T\ • 96 

7^(77 — 1) — r] 

Since we have assume that y = y c , we may use these two expressions for y to relate 77, 7r and 7 
to find that 

7T 

tt(1 - 7~ 2 ) + 7" 

This is exactly the curve along which the G 2 part of Theorem 11.31 applies . 

Finally, we may use Lemma 14.31 to invert our expression for N in terms of t. Asymptotically 



V = —i , .._ 2 - (97) 



This allows us then to express 

P(/»t(^) > ((1 " 2?7)2/ + 277(1 - rj))t - V(l - t - »7)* 1/2 a0 = P(^2(JV, M) < t), (99) 



where t is as above. It then follows after a little algebra that Theorem 11.31 applies and gives that 
these probabilities asymptotically equal 



F G 



x 



F g 



X 



(100) 



as desired to prove this part of Theorem 11.11 

Fq case: Similarly to the previous case, we rehash the height function event in terms of the LPP 
with two-sided boundary conditions event and show that the desired asymptotic probabilities arise 
from Theorem 11.31 The region of p~,p+ for which we wish to prove Fq fluctuations corresponds 
to 7T + rj > 1. We wish to prove Fq fluctuations for all y such that 1 — 2p_ < y < 1 — 2p + . 
Translating this region into rj, tt, 7 variables exactly corresponds to the region in which Theorem 
11.31 implies -Fo fluctuations. 

We wish to compute the asymptotic formula for 



P (h t (yt) > \{y 2 + l)t - 2-V3(i _ yif/h^ 



(101) 



where we have used the fact that h(y) = (y 2 + l)/2 in the region of which we are 

considering. Without loss of generality let us assume that y > (the other case follows similarly) . 
As before, set 

N — M = yt, (102) 



N + M = \(y 2 + l)t - 2- 1 /3( 1 _ y2)2 /3 fl_l 

2 v 1 + y 



The height function event we are considering has the same probability as 

P(L 2 (N,M) < t). 
Using the equations for N — M and JV + Mwe can solve for 



N 



{l + yft-l-^il-y'fl'H 1 '^ 



2\2/3*l/3, 



M = 1(1 - yft - 2" 4 / 3 (l - y 2 ) 2 /3 t V3 



X. 



These expression allow us to express 7 asymptotically as 
2 



7 



1 + 2/ 



1 + 2 2 / 3 xt~ 2 / 3 (l 



y 



2\2/3 



1 



+ 



1 



l + y) 2 (1-2/) 



+ o(f 



-2/3> 



(103) 
(104) 

(105) 
(106) 

(107) 



We may use Lemma 14.31 to invert our expression for M in terms of t. Asymptotically 

/ - + 2 1; °— — — -xM"--> + o(.V •'). (108) 



^L_ + 21 / 3 (1±^! iM v3 +( , ( mV3). 

(1 - y y (1 - y) 2 



This implies that t 2 ^ 3 



4-2/3 
(l-»)-Vs 



M 2 / 3 + o(M 2 / 3 ). This can be plugged into our expression 



for 7 and gives a new expression for 7 in term of M now: 

2 



7 



1-2/ 



1 - 2 4 / 3 xM" 2 / 3 



y 



l + y 

From this we may find that 

1-7 (1-7)7 



+ o(M 



-2/3 



1 + 7 (I+7) 5 / 3 



(1+2/) 4 / 3 . 
X M' 2 / 3 + o(M" 2 / 3 ). 



(109) 



(110) 



This can then be substituted into equation (|108|) which gives 

t = (1 + 7 ~ 1 ) 2 M + ^ + 7 ) 4/S M x / 3 x + o(M 1 / 3 ). (Ill) 

7 

Plugging this into P(L2(N, M) < t) we find that our height function probability is asymptot- 
ically equal to 

P \L 2 (N, M) < (1 + 7- x )M + iH-^M 1 / 3 ^ + o(M 1 / 3 )^ . (112) 

As it was already noted the which we are considered maps exactly onto the range of 

rj, 7r,7 which the LPP probability above is asymptotically equal to Fq(x) and hence the same 
holds for the the height function probability. □ 
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